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Rank of Projection-Algebraic Representations of Some 
Differential Operators 

Abstract. The Lie- algebraic method approximates differential operators 
that are formal polynomials of {l,x, J^} by linear operators acting on a 
finite dimensional space of polynomials. In this paper we prove that the 
rank of the n-dimensional representation of the operator 

K _ d k d k+1 d k+p 

is n— k and conclude that the Lie- algebraic reductions of differential equati- 
ons allow to approximate only some of solutions of the differential equation 
K[u] = f. We show how to circumvent this obstacle when solving boundary 
value problems by making an appropriate change of variables. We generali- 
ze our results to the case of several dimensions and illustrate them with 
numerical tests. 

1 Introduction 

The Lie-algebraic method of discrete approximations was proposed by 
Calogero in his seminal work [5], with the purpose of solving eigenvalue 
problems for ordinary and partial differential operators. Several years later, 
Mytropolski, Prykarpatsky, and Samoylenko have established a general 
algebraic-projection framework for the method that has broadened its 
applicability to a wide class of problems of Mathematical Physics [6]. The 
fundamental idea of the method is to project functions onto the space 
P n [x] of polynomials of fixed degree n and to approximate differenti- 
al operators by linear maps acting on P n [x]. In early 2000s, Bihun and 
Lustyk estimated the error of such approximation and applied the method 
to solving boundary value problems for elliptic partial differential equati- 
ons [H [2, E]- An overview of the literature on the Lie-algebraic method of 
discrete approximations can be found in pf]. 
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In this paper we prove that the rank of the n-dimensional representation 
of the operator 

K _ d*_ d k+1 d k+p 

ak dx k ak+1 dx k+1 ak+p dx k+p 

is n— k and conclude that the Lie- algebraic reductions of differential equati- 
ons allow to approximate only some of solutions of the differential equation 
K[u] = /. We show how to circumvent this obstacle when solving boundary 
value problems by making an appropriate change of variables. We generali- 
ze our results to the case of several dimensions and illustrate them with 
numerical tests. 



2 The Lie- Algebraic Method 
of Discrete Approximations 

Every differentiable function u : [a, b] — > R can be approximated by its 
Lagrange polynomial of degree n. To construct the Lagrange polynomi- 
al L[u] of a differentiable function u : [a, b] — > R, the interval [a, b] is 
partitioned by points a = xq < X\ < ... < x n = b and the values 
L[u]{xi) = u(xi) are prescribed, where i = 0, 1, 2, . . . , n. The Lagrange 
polynomial of u is uniquely defined by its n + 1 values at the points of the 
partition and can be written in the form 

n 

L[u](x) = y^ j u(x l )h(x), 

i=0 

where n 

Yi k = a i x ~ x k) 

l.{x) = 

]^[fe=o (X{ — Xk) 

is the i-th standard Lagrange polynomial, < i < n. 

Let A be a differential operator acting on differentiable functions defi- 
ned on [a,b\. We assume that A belongs to the enveloping algebra of 
the Heisenberg-Weyl algebra of differential operators {l,x, J^}; that is, 
A is a formal polynomial of {l,x, 4-}. The value of A at a function / is 
denoted by A[f]. For a differentiable function u : [a, b] — > R and a parti- 
tion a = {xo, xi, . . . , x n } of [a, 6], let u a = (w(xo), u(xi), . . . , u(x n )) T be 
the vector of the values of u at the nodes of the partition. The finite di- 
mensional representation A a : R n+1 — > R n+1 of A can be constructed from 
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the condition 

A a u a = A[L[u\] a . 

Following this condition, the finite dimensional representation Z = 
{Zkj}kj=Q of the operator d/dx can be found from the relation 



l dx 



xo) \ 



= z 



I k{xo) \ 
k{xi) 

V k{x n ) I 



where < k < n. Because lk(xj) = $kj for all k,j G {0,1,..., n}, the 
right hand side of the last equality is the k-th column of Z. Therefore, 



ry dl k 

Z jk = -( Xj ) = 



En 
m= 



m=0,m^j (xj-x m ) 



ifk = j, 



(1) 



7T& Xj Xfc 



where 7T k = YY!n=o, m ^k( Xk ~ Xm ^ for every ^ ^ {1, 2, . . . , n}. 

Because the derivative 
obtain the relation 

O) = J2 Z jk{l{x) 

i=o 

The finite dimensional representation X of the operator x of multipli- 
cation by the variable x is the square matrix of dimension (n+ 1) with the 
diagonal (xq, Xi, . . . , x n ): 



^ coincides with its Lagrange polynomial, we 

dh 
dx 



(2) 



X = Diag(.T ,.Ti, . . . ,x n ). 



(3) 



The finite dimensional representation of the identity operator 1 is the (n + 
l)-dimensional identity matrix denoted by I n+ \. 

Having constructed the finite dimensional representations of the operators 
1, x, d/dx, we can approximate every operator A in the enveloping algebra 
of {1, x, d/dx}. 



Example 2.1. Consider the boundary value problem 

f w " + w = o, 

\ u(0) = 2,u(tt/2) = 1 



(4) 



for a C 2 function u : [0, 7r/2] — > R. Let a = {xi}^ Q be a partition of 
[0,7r/2] and let Z and I n+ \ be the finite dimensional representations of 
^ and 1, respectively, for this partition. We approximate the differential 



operator A = d 2 /dx 2 + 1 by Z 2 + I n+ \ and reduce problem (0) to the 
problem of finding a vector u a G W l+1 that satisfies 

r (z 2 + / n+ i)M ff = o, , . 

\ Ho = 2, H n = 1. u 

A natural question to ask here is this: does problem (f5j) have a solution? 
More generally, when does the finite dimensional approximation of the 
original problem have a unique solution in We address this question 

in the following sections. 



3 Rank of the Finite Dimensional Representation of 
the Differential Operator J2™ =k ctp-^ 

Let 1 < k < m < n and let a&, ak+i, . . . , a m be real numbers such that at 
and a m are not zeroes. In this section we prove that the finite dimensional 
representation 

a>kZ k + cik+iZ k+1 + . . . + a m Z m 
of the differential operator 

d k d k+1 d m 

a k-n: + a k+ 1 7 U _i_l H" • • • H" Q>TTl 



' dx k dx k+l dx m 

has rank n + 1 — k. 

Lemma 3.1. Let a = {xo, x±, . . . , x n } be a partition of the interval [a, b] 
and let Z, given by formula (p]) ; be the finite dimensional representation 
of d/dx with respect to this partition. Then rankZ = n and Z n+1 = 0. 

Proof. Let us recall that the matrix Z is (n + l)-dimensional. 
First, let us show that Z n+1 = 0. Using formula p|, it is easy to verify 
that 

d k l- n 

= /\ Zk ]mjlm(x) 

for every k G N and j G {0, 1, . . . , nj. Because dxn+ { = 0, we obtain 



y ^\z n+i \ m ji m (x) — o, 



m—0 



for all x G [0, 7r/2], which, after evaluation x = Xk for /c G {0, 1, . . . , }, 
implies that Z n+1 = 0. 



Let us now prove that rank Z = n. 

There exists a nonsingular matrix B = {£>^}™ J=0 of dimension n + 1 
such that 



n 

P 



for every integer j between and n. Then 

dl n 

-^-(x) = ^ ^ Z m jl m ( K x) 

m=0 
n n 

m=0 p=0 
n 

jj=0 

Because the polynomial ^r{x) has degree n — 1, the last row of the 
matrix BZ must be zero. Therefore, rankZ = rank£>Z < n. 

It remains to show that rank Z cannot be strictly smaller than n. 
Suppose, on the contrary, that rank Z < n. Then there exist n linearly 
dependent columns in Z . Without loss of generality we assume that the 
first n columns of Z are linearly dependent; denote these columns by 
t4q, Ai, . . . A n _\. There exist constants cko, ai, . . . , a n -i, not all zeros, such 
that 

n— 1 
3=0 



But then 



n— 1 ^ n— 1 n 

j=0 j=0 m=0 

n n—l 

771=0 7=0 

= o, 

which, after integration, gives us 

77-1 

3=0 
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where C is a constant. 

By setting x = Xj in the previous equality, where j G {0, 1, . . . , n — 1}, 
we obtain a® = a\ = . . . = a n -i = C ^ 0. On the other hand, if we set 
the very same equality, we obtain C = 0, a contradiction. We 
have thus proved that rankZ cannot be smaller than n. Q.E.D. 

Lemma 3.2. Let H be a nilpotent square matrix of dimension n + 1 such 
that ranki7 = n and H n+1 = 0. Then ranki7 fc = n + 1 — k for all 
k G {0, 1, . . . , n}. Moreover, the Jordan form of H is 





where I n is the identity matrix of dimension n. 

Proof. Let J be the Jordan form of H so that H = B l JB for some 
invertible matrix B and diag(J) = (Ao, Ai, . . . , A n ), where A^ are the ei- 
genvalues of H. 

Then = H n+1 = B^J^B and J n+1 = 0. Because diag(J n+1 ) = 
(A£ +1 , Af +1 , . . . , A; i+1 ), the matrix J has only one eigenvalue A = of 
multiplicity n + 1 and is given by formula ([6]) . 

It is now easy to verify that 

I n +i-k 



J = U o J « 

and rank J k = n+l — k. Thus, rank H k = rank J k = n+ 1 — k as required. 
Q.E.D. 

Theorem 3.3. Let P(z) = akZ k + ak+iz k+1 + . . . + a m z m be a polynomial 
of degree m with real coefficients, where k G {0, 1, ... , m} and a k ^ 0. // 
a square matrix B is nilpotent, then rankP(5) = ranki^. 



Proof. 



Let us rewrite P(B) as P(B) = a k B k (I + ^B + ^B 2 

a^ B m-k\ De £ ne £, + ^+1 #2 + _ _ _ + a^gm-k^ 

a k ' a,k a>k a k 

Because the matrix B is nilpotent, D is nilpotent. But then the matrix 
I+D is invertible. Indeed, the series Y^=o(~ D) p is a finite sum and equals 
(I + D)~\ 

In summary, P(B) = B k A, where A is an invertible matrix. Therefore, 
rankP(5) = rankfi fc as required. Q.E.D. 



The next theorem follows immediately from theorem [3^3] and lemmas |3T 
and [3T2 
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Theorem 3.4. Let o = {xo, Xi, . . . , x n } be a partition of the interval [a, 6] 
and let Z , given by formula be the finite dimensional representations 
of the operators d/dx and x, respectively, corresponding to this partition. 
Let P(z) be a polynomial defined in the hypothesis of theorem 3.3. Then 
rankZ = n, Z n+1 = 0, and rankP(Z) = rankZ^ = n + 1 — k. 



It is natural to ask whether theorem |3.3| can be generalized to the case 
where the coefficients a p , p G {k, k + 1, . . . , m}. of the polynomial P(z) = 
cikZ k + ak+\z k+l + . . .+a m z m 1 are polynomials in x. In other words, suppose 
that a p , p G {k, k + 1, . . . , m}. are polynomials, matrix B is nilpotent and 
matrix X is diagonal. Is it true that the matrix P(X, B) has the rank of 
B k l The next example shows that, in general, this is not true. 

Example 3.5. Let 



B= . „ , and X = 





It is easy to verify that B 2 = 0, hence B is nilpotent. Consider the 
polynomial P(x, z) = 1 + xz. The matrix P(X, B) = / 2 + X.B does 
not have the full rank if b — a = — |. Indeed, 



1 — 2a —a 
46 1 + 26 



/ 2 + = 

and det(/ 2 + XB) = 1 + 2(6 - a) = if 6 - a = 



Let us revise example |2.1| in the view of what we have proved. The 
finite-dimensional approximation of the differential equation u" + u = 
on the interval [0, 7r/2] is given by 

Z 2 u a + u a = 0, (8) 

where u a G WL n is a vector whose coordinates approximate the values of 
the solution u at the nodes of a partition a of [0, 7r/2]. 

By the last theorem, the rank of the matrix Z 2 + 7 is n + 1, thus the 
system of linear equations ([8]) has a unique solution. On the other hand, 
we know that the space of solutions to the equation u" + u = is two- 
dimensional. Therefore, such direct application of the Lie-algebraic method 
does not help us approximate all the solutions of the differential equation, 
only some of them. It is possible, however, to still approximate solutions of 
initial and boundary value problems using the Lie-algebraic method. We 
simply need to make an appropriate substitution so that to ensure that 
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the solution of the differential equation automatically satisfies initial or 
boundary conditions. 

For example, to solve the boundary value problem 

f u" + u = 0, x G [0, f ] . , 

\«(0) = 2,«(f) = l. iyj 

we make a substitution u = (2 — 2/irx) (x(x — ir/2)v + l) and obtain the 
equation p(x)v" + q(x)v' + r(x) , u = s(x), where 

2 

= X 3 + 3x 2 — 7TX, 

7T 

6 

q(x) = 2( x 2 + 6x — 7r), 

7T 

/ \ 12 „ /„ 2 7T 

r(x) = x + 6 + x(2 x)(x— — ), 

7T 7T 2 

2 

s(x) = —x — 2, 

7T 

whose discrete approximation is 

(p(X)Z 2 + q{X)Z + r(X))v ff = s(X), 

where a is a partition of [0,7r/2]. 

We have chosen the partition of the interval [0.001, |] with equally 
distributed n nodes X{ = 0. 00 1 + — 0.001). We compared the values of 
the vector v a with the values (u(xq),u(xi), . . . u(x n )) of the exact solution 
u(x) = sinx + 2cosx to the boundary value problem ([9]). The values of 
the errors 

n 

and 

E max (n) = max I - u(x k ) I 

fce{0,l,...,n} 

are given in table 1. 

We have compared the results obtained by means of the Lie-algebraic 
method with the results obtained using the standard shooting method. 
Following the shooting technique, we have split boundary value problem (|9]) 
into two initial value problems: 

{ w " = ~ W (10) 
\ w(0) = 2,u/(0) = 0, W 
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and 

r v" = -v 

{ v(0) = 0,t/(0) = 1. 
The solution of problem (|9j) can be found by the formula 



it = w + 



1 



7T 



«(!) 



■U. 



We solved the initial value problems (10) and (11) using finite di- 
fferences. 

The error of the shooting method for problem ^ is given in table 1, 
which demonstrates fast convergence of the Lie-algebraic scheme. 



Method 


n 


4 


8 


12 


16 


Lie-Alg. 


E(n) 


2.2788e-04 


9.5522e-07 


3.9033e-09 


1.5205e-ll 


Lie-Alg. 




1.1466e-04 


2.5575e-07 


6.8542e-10 


1.9955e-12 


Shooting 


E(n) 


2.71e-02 


1.39e-02 


9.3e-03 


7.0e-03 


Shooting 




l.le-02 


2.7e-03 


1.2e-03 


6.7013e-04 



Table 1. The errors E(n) and E max (n) of the Lie-algebraic scheme and the shooting method 
for boundary value problem Q, where n is the number of subintervals in the partition of 
[0.001, 5] for the Lie-algebraic scheme and [0, \\ for the shooting method. 



4 The Lie-algebraic Approximations 
in Several Dimensions 

4.1 Lagrange Polynomials of Several Variables 

In this section we will generalize the results proved in the previous section 
for the case of several dimensions. More precisely we will study the exi- 
stence and the number of solutions of the Lie-algebraic discrete approxi- 
mations of partial differential equations. 
Consider a differential equation 

Au = f, 

where Q is a (i-dimensional cube [a±, b\} x [a 2 , ^2] X . . . X [a^, bd] and / : 
Q — > R belongs to the domain of the differential operator A. We denote 

,'Y> ( ry 1 ryd\ 

tXj I tX/ • ■ ■ ■ ■ >/ / • 

We partition each interval [a^, where 1 < k < d, with the points 
Ofc = x k o < x k i < . . . < x h nk = bk, and denote the basic Lagrange 
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polynomials with respect to this interval by 

n k n k 

i k j(x k )= n (x k ~x k m )/ n (^-x fc m ), 

m=0,m^j m=0,m^j 

where j = 0,1,.. . ,n k . 

The basic Lagrange polynomials on Q are defined by 

I J (i 1 ,i 2 ,—,id)( X i X , • • • , x ) = ^ h{ x )l ijy x ) ■ ■ - l i d { x )■ 

For brevity we introduce the multi-index notation (i) = {i\, i<i, • • • id), 
where < ij~ < n k for every k = 1,2, ...,d. The set of all multi- 
indices is denoted by 3 and the nodes of the partition of Q are x^ := 
(x 1 ^, x 2 i 2 , . . . , x d i d ), where (i) G J. In this notation, the basic Lagrange 
polynomials on Q can be written as L^(x), where (i) G J, and the 
Lagrange polynomial of a function u : Q — > E is 

(i)ea 

Let n a -dimensional square matrix denote the representation of the 
differential operator ^ with respect to the partition {x a k}^l of [a a , b a ], 
where a = 1, 2, . . . , d. It is easy to verify that for 1 < a < d 

d -L {i) (x) = Y,W ia \ Mi) L {j) (x), 



where 



dx 



\ji,j2,-,jd),{ii,i2,-,id) ^ f /. fc«o' 



£ is the Kronecker symbol, and a "hat" above the index i a means that the 
index i a is omitted. 

We consider as a linear map from V := M ni+1 ®M n2+1 ® . . 
to itself, which acts in the following way. Let {e k j}™= be the standard basis 
of R nk+1 . Then 

W^\e\ ® e\ ® . . . ® e d ld ) = ^ ^, {j /h ® ^ ® ■ ■ ■ ® e V 

We identify the map W^ a ) with a square matrix of dimension 
7V:= (m + l)-(n 2 + l)-...-(n rf + l) 
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that represents the linear map in the standard basis of HH N = V as 

follows. We enumerate the multi-indices in 3 by introducing the following 
bijection * : J — >• {1, 2, . . . , N}: 



*(o,o,o,. 


..,o) 


= 1, 


*(i,o,o,. 


-.,0) 


= 2, 


*\Jll, U, U, . 


. . ,UJ 


n 1 1 

— Til + J-, 


*(o,i,o,. 


..,0) 


= (m + i) + i, 


*(m,i,o,. 


..,o) 


= (m + 1) + (m + 1), 


*Oi,n 2 ,0, . 


..,o) 


= n 2 -(ni + l) + (m + l), 


*(0,0,1,. 


..,o) 


= l-(n 2 + l)(ni + l) + l, 


*(ni,n 2 ,n 3 , . . 


• ,n d ) 


= n d ■ (n d _i + 1) . . . (ni + 1) 



+ . . . + n 3 ■ (n 2 + l)(m + 1) + n 2 ■ (n 2 + 1) + m + 1. 

It is easy to see that 

= *(ii,i 2 , . . . ,i d ) 

= i d (m + l)(n 2 + 1) . . . (n d _i + 1) + . . . + z 2 (ni + 1) + zi + 1 
for all (i) e 3. 

Let {s'i}^! be the standard basis of R^. For every multi-index (i) G 3, 
we identify the basis element with the basis element e 1 ^ <S> e 2 i 2 <S> ... <S> 
e d i d of V. Using this identification, it is now easy to compute the matrix 
that represents the linear map : V — > V: its *(i)-th column 

equals 



W {a) g Hl) = W^ a \e l h ®e 2 i2 ®...®e d id ) 

M) 



Therefore, 

*W>*0) W.O) h,...,i a ,...,i d %c " 3a v 7 

In the following, we identify iy( a ) and 

11 



Proposition 4.1. The matrix representation of the differential operator 
-J^ is given by the N x N -dimensional matrix 

W {a) = I ni+1 <8> . . . ® I na _ 1+1 <g) (2) (8) ... (8) (13) 

where <S> denotes tensor product and I nk +i is (nk + 1) -dimensional identity 
matrix. 

We prove the above proposition for the case d = 2 in the next subsecti- 

d 

on, where we discuss the structure of matrix representations of — — for 

dx a 

the two-dimensional case. The proof for multi-dimensional case (d > 2) is 
similar and is left to the reader. 

It can also be shown that the finite dimensional representation of the 
operator of multiplication by a variable x a equals 

= / ni+ i <8) . . . <8> J„ a _ 1+ i <8> X (a) I na+1+1 (g) ... (g) J n<J+1 , (14) 

where X^ a ^ = Diag(x a o, x a i, ■ ■ ■ , £ a n a )- 

We end this subsection with stating useful identities for matrix W^ a \ 

Lemma 4.2. The matrix 

[W (a) ] k = In x+ 1 ® ■ ■ ■ ® ® t^ (a) ]" ® 4 a+1+ l ® • ■ ■ ® (15) 

represents differential operator g^k for all k G {1,2,..., n a } ; and 

lyW.iy^) = / ni+1 ^ . . . ^ / na _ i+1 ® ® / na+1+1 ® . . . (16) 

represents differential operator d ° a g x n for all a, (3 G {1, 2, . . . , <i}. 

Moreover, the matrices W^ a \ where a G {1, 2, . . . , d}, commute and 
are nilpotent. 



Formulas (15) and (16) follow from a composition property of tensor 
products: For every pair A, B of linear maps (between linear spaces) , 
(A <g> Bf = (A (g) B)(A ® B) = A 2 ® 5 2 . The last statement of the 
last lemma is an immediate consequence of equations (15) and (12), and 
lemma 13.1 
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4.2 Matrix Representations of ^ 
for the Two-dimensional Case 

Let us observe the structure of the matrices and in the two 

dimensional case (d = 2): 

( ... \ 

... 



\ 







zw 



J 



(17) 



where the number of blocks Z^> on the diagonal of is n<i + 1; 



7(2)7- 



' ^0,0^1+1 

7(2) r 7(2) r 

Z 'l,0 i ^i+1 Zj l,l 1 n 1 +l 



7 (2) /" 



\ 



, A(2) . A(2) . 



0,n 2 J ni+l 
Z 'l,n 2 i n 1 +1 

7 (2) r 
^n 2 ,n 2 J ni+l / 



where is the (ni + l)-dimensional identity matrix. 

Let us recall that for two linear operators A : X — > X and B :Y — > Y 
acting on linear spaces X and Y, their tensor product A <g> B : X <g> Y — > 
X <g> Y is defined by A <g> B(x <g> y) := Ax <g> By. 

Let {eo, . . . , e Hl } and {fo, . . . , f n2 } be the standard bases of and 
M n2+1 , respectively. Let N = (ni + l)(n2 + l) and (i) = (21,22) £ 3. Having 
identified the *(z)-th vector in the standard basis {gi,g2, ■ • • , 9n} of 
with e ix (g) /j 2 , it is now easy to observe that = Z^> (g) 7 n2 +i and 

W {2) = /„ 1+ i ® ^ (2) - 

Indeed, the *(z)-th column of the matrix TJ^ I n2 equals 

[Z {1) ® I n2+l ]. Al) = ®I n2+1 {e h ® f i2 ) 

= Z^e H ®I n2+1 f l2 



ii=o 

ni n 2 



zi e j'i 



'2 



h =0.72=0 

= E ^i(i),*(i)^0> 
(i)eJ 
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.72 



which is the *(i)-th column of The equality = I Ul +i <8 can 
be verified in a similar fashion. 

4.3 Rank of Matrix Representations of Operators ^ 
and Their Compositions 

In this subsection we compute rank of matrix representations of operators 

l3l 



d 

x 

We start with proving the following simple fact. 



, , , and their compositions, by generalizing results of section 



Lemma 4.3. Let A : E TO — > E m and B : W — > W be linear maps with 
ranks mi and respectively. Then rank(^4 <g> B) = mi • m-2. 

Proof. Denote 

kA '■= dim(ker(74)) = m — m\ and k B '■= dim(ker(5)) = p — m<i- 
Let {f\, f2, ■ ■ ■ , fk A } be a basis of kei(A); complete it to a basis 

{/l? fli • ■ ■ ? fkAi fkA+li • • • 1 fm} 

of M. m . Similarly let {gi,g2, ■ ■ ■ , gk B } be a basis of ker(5); complete it to a 
basis {#1, g 2 , . . . , g kB ,gk B +h ■ ■ ■ , 9 P } of M p . For every u G M. m and i> G E p , 
A (8 B{u (8 v) = Au <S> Bv = if and only if Au = or Bv = 0. Therefore, 



ker(A^B) = span{/i <8> / 2 <8> g%, ■ ■ ■ , /k A ® & : i = 1, 2, . . . , k B } 

+ span{/i <8> / 2 <8> . . . , f kA <8> g { : i = k B + 1, . . . 7 p} 
+ span{/ z <8 #1, /* <8> #2, • • ■ , fi <8> #fc B : i = &a + 1, . . . , m}, 

and 

dim ( ker(^4 (8 5)) = &a • &b + ■ ^2 + k B ■ m\ = pm — mi?7i2, 

which implies rank(.4 B) = m\m2 as required. Q.E.D. 

Using the previous lemma, it is easy to verify the statement of the 
following proposition. 

Proposition 4.4. Let B and C be (n a + 1)- and (np + 1)- dimensional 
square matrices, respectively. Let 

M := I ni+1 (g) ... (8) I na _ l+ i (8) B (8 <8> . . . <8> 

and 

K := 4 1+ i(8>. • .<8>/ nQ _i+l®£®Jn„ f i+l®- • ■®4«_i+l®C^8>/», 8+1 +i<8>. • -®4 d +l- 
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Then 



and 



rank M k = rank B k 



fni + l)...(n d + r 



rank A' = rank B ■ rank C 



(n a + 1) 

;m + 1) . . . (n d + 1; 

(re a + l)(re^ + l) 



From the previous proposition, formulas (13), (15), (16), and lemmas 3.1, 3.2 
we conclude the following. 

Proposition 4.5. Let a, f3 G {1, 2, . . . , d}. For all 1 < k < n a and 

1 < / < fi/3, the rank of the matrix representation [W^] k of 

;rei + l)(re 2 + l)...(re d +l) 



d(x a ) k lS 



rank[W (a) ] fc = (n a + l-k)- 



[n a + r 



an 



Q(k+l) 



IS 



d the rank of the matrix representation \W^] k \W^] 1 of g^^gn^y 

rank[^>]^<«]< = (^ + 1 - fc) (^ + 1 - Z) ( ^ + D(n, + 1) . (», + 1) 

(n a + l){np + l) 

Moreover, the matrix is nilpotent: [H/( a )] n « +1 = 0. 

Let a = {x^}i e j be a partition of Q, = [ai, &i] X . . . X [a d , 6 d ] as before. 
Denote the number of nodes in the partition by N. Recall that the matri- 
ces and denote the matrix representations of the differential 
operator and the operator of multiplication by x a : respectively, with 



respect to the partition a (see equations (13) and (14)). Consider partial 
differential equation 



P(x , . . . , x , 1 , . . . , j)u /, 



(19) 



'' Ms a formal polynomial of its arguments. 



dx 1 ' ' dx d 

where P(x 1 , . . . , x d , Jr, . . . , ^) 
Let F = (/(x( 7l) ), /(x (72 )), . . . , /(x (7jv ))) 5 where 7fc G J and *(7 fc ) = k for 
all G {1,2,..., N}. 



To solve PDE (19) using Calogero method, we need to solve the system 
linear equations 

P{X^\ . . . , X®, W {1 \ . . . , W {d) )U = F (20) 

for U G WL N ; the i-th component Ui of the solution vector U, if it exists, 
approximates the value u(x^), where #(7) = i, of solution w of equati- 
on (19) at the node Xhy\ of the partition a. It is therefore important to 



know the conditions under which the system of linear equations (20) has 
a unique solution. 
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Theorem 4.6. Let P(z 1 , z 2 , . . . , z d ) be a polynomial with real coeffici- 
ents. The matrix P(W^\ W^ 2 \ . . . , W^ d >) has full rank if and only if the 
constant term of P is nonzero. 

Proof, (i) If the constant term of P, which we denote by a, is 
nonzero, then P(z 1 , z 2 , . . . , z d ) = a + Pi(V, z 2 , . . . , z d \ where P\ is a 
polynomial in (zi, . . . , zj) whose constant term vanishes. Because matri- 
ces W^ a \ where a G {1,2, commute and are nilpotent (see 
lemma 4.2), the matrix Pi(W^\ . . . , W^ d >) is nilpotent. But then the 



matrix I N + \Pi{W^ l \ . . . , W^) is nonsingular, as the formal series 
expansion of (In + ^P\(W^ l \ . . . , M^^)) _1 is a finite sum. Therefore, the 
matrix 

P(W {1 \ . . . , W {d) ) = a(I N + -Pi(W {1 \ . . • , W {d) )) 

a 

has full rank as required. 

On the other hand, if the constant term of P vanishes, then the matrix 
P(W( 1 \ . . . , W^) is nilpotent, and therefore does not have full rank. 
Q.E.D. 

Example 4.7. We finish this article with an application of the Lie- 
algebraic method to the following partial differential equation of hyperbolic 
type [7] 

u xx - u yy + yu x = 4(y 2 - x 2 ) sin(l - x 2 - y 2 ) - 2xy cos(l - x 2 - y 2 ) (21) 

on region Q = {(x, y) : x 2 + y 2 < 1} with the boundary condition 

u\ d n = 0. (22) 

This problem has solution u(x,y) = sin(l — x 2 — y 2 ). We compute the 
approximate solution of the problem using Lie-algebraic approximations 
and compare our results with the exact solution. 

To incorporate boundary conditions, we make the change of variables 

u = (1 — x 2 — y 2 )v 

, which yields the partial differential equation 

(1 - x 2 - y 2 )(v xx - v yy + yv x ) - 4xv x + 4yv y - 2xyv (23) 
= 4(y 2 — x 2 ) sin(l — x 2 — y 2 ) — 2x?/cos(l — x 2 — y 2 ). 

We partition the square D = [—1,1] x [—1,1] that contains region Q 
with points a = {(xi,yj) : < i < ni,0 < j < 77,2}, where o x = {^}^= 
and G y = {yi}™=Q are partitions of the interval [—1,1]. 
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Let Z x and Z y be matrix representations of d/dx on [—1, 1] with respect 
to partitions o x and o yi respectively. Similarly, let X and Y be matrix 
representations of the operator of multiplication by x on [—1, 1] with 
respect to partitions o x and o yi respectively. Matrices Z x and X have 
dimension n\ + 1, while matrices Z y and Y have dimension + 1. Finally 
the identity matrices / ni +i and / n2 +i represent the identity operator with 
respect to partitions o x and a y: respectively. 

The following table gives matrix representations of the operators d/dx 7 
d/dy, x, y, and 1, with respect to partition a. 



Operator 


d 
dx 


d 

dy 


X 


y 


1 


Repres. 


Z x <S> In 2 +l 


Irn+l ® Z y 


X <g> I n2 + 1 




Inx+l ® In 2 +1 


Notation 


Z x 


A 

Zy 


X 


Y 


I 



Table 2. Finite dimensional representations of basic differential operators in two dimensions. 



Denote the right hand side of equation (23) by 

y 2 ) — 2xycos(l — x 2 — y 2 ). 



f(x,y) = 4(y 2 - x 2 ) sin(l 



x 



Equation (23) is represented by the system of linear equations 

Kv a = F, 

where 

K = {I -X 2 - Y 2 ){Z 2 X - Z 2 + YZ X ) - 4XZ X + 4YZ y - 2XY 



(24) 



(25) 



F = (f(x , y ), f(xi, 2/0 ), • • • , f(x nv y ), f(x , y{), f(x u y{), . . . , f(x Hl , y n2 )) J \ 
and the unknown vector v a E R( n i +1 )( n 2+i)_ 



After solving system of linear equations (24), we compute the vector u c 
of approximate values of u as follows: 

u a = (I - X 2 - Y 2 )v a . 



We have compared the approximate solutions with the exact solution by 
computing the maximal error E{n\,ri2) = maxvywg l^o-Wij) — u(xi,yj)\ 
and the average error E a {n ll n 2 ) = {ni+l ] [n2+l) E(i,j)gj \W]*{i,j)-u(x h yj)\. 
The values of E and i? a for different number of partition points are given 
in table 3, and the plot of the approximate solution u a is given in figure 1. 
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(ni,n 2 ) 


(10,10) 


(15,15) 


E(ni,n 2 ) 


0.0064 


0.002 


E a {n u n 2 ) 


2.56e-04 


2.63e-05 



Table 3. Maximal and average error of the Lie-algebraic method for problem (21), (22) 




Figure 1. Plot of the approximate solution u a to problem (21), (22) 
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